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I. INTRODUCTION 

The second quantization is the effective method, which is used for description of physical processes in quantum 
theory of many-particle systems. The key role of this method consists in introduction of creation and annihilation 
operators of particles in a certain quantum state. The operators of physical quantities are constructed in terms of 
creation and annihilation operators. Such a description of quantum many-particle systems implies the particles to 
be elementary (not consisting of other particles). Moreover, it is absolutely accurate despite of the possible existence 
of compound particles. Since the interactions between particles may lead to formation of bound states, the standard 
second quantization method becomes too cumbersome. For this reason the construction of approximate quantum- 
mechanical theory for many-particle systems consisting of elementary particles and their bound states represents an 
actual problem. In this theory it is necessary to introduce the creation and annihilation operators of bound states as the 
operators of elementary objects (not compound). In addition we should preserve the required information concerning 
internal degrees of freedom for bound states. Such an approach may be actually realized in the approximation in 
which the binding energy of compound particles is much greater in comparison to their kinetic energy. 

The problems leading to the necessity of accounting for compound particles along with elementary ones are typical 
in studying of the interaction of irradiation with matter (the matter consists of neutral atoms or molecules). The 
atoms or molecules may be both in the ground and excited states when it is necessary to take into account their 
internal structure retaining the convenience and simplicity of the second quantization method. Such a situation 
should take place in describing the experiments on laser cooling of atoms 1-3 or experiments on observation of Bosc- 
Einstcin condensates in magnetic traps, where a condensate temperature is reached by laser cooling of gases of alkali 
atoms 4 . The similar problem of accounting for the bound states occurs also in quantum chemistry for example, in 
studying of chemical reactions. 

Thus, the construction of the quantum theory of many-particle systems in the presence of bound states of particles 
within the second quantization method represents not only academic but also applied interest. The formulation of the 
quantum theory itself is to be done so that the operators of physical quantities are constructed through the creation 
and annihilation operators of elementary particles and their bound states. 

In this paper we demonstrate a possibility of construction of the above-mentioned quantum-mechanical theory by 
considering a system of two kinds of charged fermions. In this theory the compound particles (atoms) consist of two 
fermions of various kinds. The choice of such model is not associated with principle difficulties and makes it possible 
to simplify calculations and to obtain the visual results. 

As applications of the developed theory we study a spontaneous radiation of two-component excited atom and 
obtain the expression for its probability (see, e.g., 5 ). We also investigate the attractive forces acting between neutral 
atoms at long distances (the van der Waals forces) and derive the expression for the potential of the van der Waals 
forces 6 . Finally we consider a question concerning the scattering of photons and elementary particles by bound states . 
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II. CREATION AND ANNIHILATION OPERATORS OF BOUND STATES OF PARTICLES 

Consider a system consisting of two kinds of fcrmions with masses mi and m 2 . As was mentioned in the Introduction 
this case is more visual in order to demonstrate the procedure for obtaining the operators of basic physical quantities 
within the second quantization method in the presence of bound states of particles. We shall see that the developed 
methodology does not contain any difficulties for its generalization (the bound states of more than two particles, the 
presence of bosons in the system, the bound states of bosons and boson- fermion bound states) . For simplicity we do 
not take into account a spin variable (this can be done without any difficulties). 

Let t/>i(x), ?/> 2 (x) be the annihilation operators of two kinds of fermions at the point x, 

V'i(x)|0) =&(x)|0) =0, 

where |0) is a vacuum state vector. Then the state vectors 

|xi, . . . ,x„,yi, . . . ,y m ) = t/4 (xi) . . .-0j f (x n )^(yi) . . . ^ (y m )|0), (1) 

(n,m = 0, 1,2, . . .) form a basis in the space of states H. In these states the particles are at the certain points 
Xi, . . . , x„; yi, . . . , y m G R of the coordinate space. The state vectors Q satisfy the orthogonality and normalization 
relations and form a complete set of state vectors. 

We assume that the particles of two different kinds form a bound state specified by the wave function 

I \ I \xf v\ v ™ixi+m 2 x 2 

<yS Q (Xl,X 2 ,x) = <£ Q (Xl - X 2 )0(X - X), X= ; , (2) 

mi + m 2 

where x is the space coordinate and a are the quantum numbers of the bound state (atom) (we suppose that the 
particles of the same kind do not form the bound states) . The corresponding state vector has the form 

|a,x) = y ebcidx 2 </? Q (xi - x 2 )<5(x - X)^(xi)^(x 2 )|0). 

For this reason the operators 

<^>+(x) = J dxidx 2 <p a (xi -x 2 )6(x - X)^" (xi)^" (x 2 ) (3) 

we shall call as the creation operators of the bound states (atoms) , so that 

0+(x)|O) = |a,x>, &,(x)|0>=0. 

If the atom has a certain momentum, then its state vector is given by 

Kp) = J dxidx 2 v3 Q (xi - x 2 )e lpX '0 1 f (xi)V^ (x 2 )|0), 

where V is the volume of the system. The corresponding creation operator 0„(p) of the atom in the state with 
momentum p is defined by 

hp) =0+(p)|O>, 0+(x) = -L^^+(p) e -^. 



Taking into account that 

J dyi<p* a (yi - y2)<Pfi(yi -y2) = S a0 , 

it is easy to find the following commutation relations: 

[^a(x),<y5+,(x')] = 6 aa '6(x.- x') +Xao'(x,X / ), [0 a (x) , <f a , (x')] =0, 

where 

W(x,x') = |d y d y V*(yWyO{^(x+^y)^(x' + ^y')5(y-y'-^(x-x' 

+4+ (x' - ^y') ^ 2 (x - ^y) S (y - y' + ^(x - x')) }, M = m 1 + rn 2 , 
moreover, x aa /(x, x')|0) = 0. The vectors 

n m I 

i X1 , . . , yi , . . , Z1) . . .) = n^ + (x 4 ) n ^ + (y fc ) n &>*)io> ^ 

^^^^^T^ 1=1 fc=1 i =1 

have an obvious physical meaning under the following conditions: 

\xi - Xj | > a, |yi - y,-| > a, \z t - Zj\ > a, 

l x * -yj\Z l x i - z jI <: a , |y« - z jI £ a > ( 5 ) 

(x, y,z € i? a , a >> ro, ro is the radius of the bound state; the definition of a see below). In this case the elementary 
particles and their bound states are at the certain space points. 

Notice that the state vectors do not form a basis in the Gilbert space H if (JSJ are valid. However, their linear 
span, which is a totality of the following vectors 

J J dyi... J dz 1 ...C(x^^,y^ (6) 



form a subspace H a of the space H . Let us show that the state vectors (@J [with conditions (JSJ)] form an orthonormalized 
basis in the subspace H a . To this end it is necessary to take into account that in calculating of the following vacuum 
averages: 

(0|&( Xl ) . ..&(x 3 ) . ..0„(x)0+ (x') ■ • .4 + (*' 2 ) • - .#(xi) ■ • ■ |0) 

we can use the Wick theorem with the following contractions: 

^(x)^+(x') = (0|^(x)7/;+(x')|0) = <W5(x-x'), ^(x)Vv(x') = °. 

if to consider the operators ifa, ip referred to the moment of time +0, and operators <p + , tp^, ip\ to the moment 
of time —0. In addition, we should remember that the creation and annihilation operators ip a , depend on V^t/v 
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(see 0). We also assume that the wave function J2J of the atom differs from zero for |xi — X2I < tq. Taking into 
account Q3J and noting that for — x 2 | > a 

^ Q (z)4 + (x 2 )^+(x' 1 ) = f dz x dz^ a (zi - z 2 )*(z - Z)^ 1 (z 1 )^ 2 (z 2 )V-+(x 2 )^+(x' 1 ) = ^(xi - x 2 ')<5(z - X') = 0, 

ab b a J abba 



, _ TO i x 1 + m 2 x 2 
mi + fn-i 



and also 



<p a (z)<p+,(z') = / dzidz 2 dz l 1 dz 2 ip* a {zi - z 2 )6(z - Z)(p a >(z' 1 - z 2 ')8{z - Z')tpi(z 1 )tp 2 (z 2 )ip 2 (z' 2 )i>t( z i) 

ab ba J abba 

= J dz 1 dz 2 if* a (z 1 - z 2 )S(z - Z)( y 9 Q /(zi - z 2 )S(z' - Z) = <S QQ '6(z - z') 
(the double contractions correspond to the operators ip a , </?+) we get 

(o| ^( Xl ) . . . Myi) (zi) • • • vt (zi) . . . ^ 2 + ( yi ) . . . MO • ■ ■ |o) = 

v >v >v y 1 ^ - /v > 

V V V v > V V 

n m / m' n' 

= ^m^mm^^^^^ ^ X l ~~ x 'l ) ■ • -Jfa ~ Yl) ■ ■ / ( Z l ~ Z 'l)<Wai • • ( 7 ) 

n m 1 

This relation shows that the vectors © form the orthonormalized basis in the subspace H a if to consider the creation 
and annihilation operators (p^(z), (p a (z) as the Bose operators, which commute with ^j(x), -0^"(x). The quantity V x ' 
in Q) is equal to +1 if the number of permutations of the arguments x' x . . .x' n is even and it is equal to —1 if the 
number of these permutations is odd. The quantity V y < is defined similarly. However, the case when it is necessary 
to take into account a more complicated arrangement of the contractions between ip a , tp+ : 

X4X4 . . . X 2 X 2 . . . Xixi . . . X3X3 . . . 

a b c d e e d c f, a 

can take place (the odd indices correspond to the creation operators and even to annihilation operators). Since 
x' a = x' 2 , x 2 = x 3 , x 4 = x 3 , 

l x i - x il ^ r o, |xi-x 2 |<r , |x 3 -x 2 |<r , |x 4 -x 3 |<r , 

whence |x4 — xi| < 4ro. Therefore, if a < 4ro, then the above contractions are to be taken into account. If however, 
we consider the n-particle states, then the discussed arrangements of the contractions are not to be taken into account 
for n < a/ro (we emphasize that rg/a is a small parameter in our problem, see pi|l). The said above concerns to 
the majority of problems in non-relativistic quantum theory, where a finite number of particles is usually studied (see 
sections V, VI). 

With the use of Q it is easy to find the projection operator Vn a onto the subspace H a : 

k-\-m-\-l<n 
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such that the operators A acting in the subspace H a £ H correspond to the operators of physical quantities A that 
act in H (hereinafter the sums of 1, 2, . . . rt-particle subspaces over the bound states are considered; n ~ a/ro >> 1), 

A = r Ha AV Ha . (8) 

Let us introduce now an auxiliary space H with the Fermi creation and annihilation operators \i ( x ) > X.t ( x ) > 
Xi(x),X2(x) and Bose creation and annihilation operators 77+ (x), ?) a (x) and take the vectors 

|xi, . . . , yi, . . . , zi, . . .) = xt ( x i) • ■ ■ xt(yi) ■ ■ ■ Vti z i) ■ ■ ■ |0), 

as a basis of this space, where |0) is a vacuum vector in H. Then the linear span of the vectors 

|xi, ...,yi,...,Zi,...) £ H a , x,y,z £ R a 

determines the subspace H a of the space H. 

Now we can easily establish the isomorphic correspondence between H a and H a : 

|xi, . . . ,yi, . . . ,zi, . . .) |xi, . . . ,yi, . . . ,zi, . . .), 

which preserves the scalar product 

(x', . . . , y', . . . , z', . . . |x, . . . , y, . . . , z . . .) = (x', z 7 , ... |x, ... ,y, ... ,z .. .), x,y, z; x',y', z' £ R a . 

We can also establish the isomorphism between the operators A •<=>■ A, acting in the spaces H a and H a according 
to the formula: 

(x\ . . . , y', . . . , z', . . . \A\x, . . . , y, . . . , z . . .) = (x', . . . , y', . . . , z', . . . |I|x, . . . , y, . . . , z . . .). (9) 

This isomorphic correspondence is remained after the multiplication of the operator by a number, after addition of 
operators, and after multiplication of operators: 

XA ^=> XA, A+B ^=> A + B, AB ^> AB. 

The formulas © , © lead to 

(xV..y,...y,...|I|x,...,y,...,z...) = x' y' ■/ 1 x y /..... (10) 

This relation determines the operators of various physical quantities A acting in H a and, hence, transfers the quantum 
theory in which the compound (the bound states) and elementary particles exist on an equal basis from the space of 
states H onto the space of states H a . We would like to remind here that entering is determined by (J3J. 

The relation (|10|l defines the operator A in H a uniquely, but it does not define it uniquely in H . It is evident that 
the operator A acting in H (continued from H a to the whole space H) is determined up to the term A' , the matrix 
elements of which are zero in the space H a (A = A' + A"). If to introduce the projection operator onto the 
subspace H a and to require that the operator A has no nonzero matrix elements in the orthogonal subspace, then 



the operator A will be defined in H uniquely, A = Vjj A"T'h . We will omit the projection operator Vfj when 
constructing the operator A acting in H. The reason for this is the assumption that the matrix elements of the 
operators (in the position space), corresponding to a quite large external parameter R ~ |x, — Xj\ give a dominant 
contribution to the quantum- mechanical processes. Further we will assume that R » r (usually in case of particle 
collisions R^ 1 ~ ym£, where £ is the particle kinetic energy). The latter inequality makes it possible to choose 
parameter a (see |[SJ) as follows 10 : 

R» a» r . (11) 

The matrix elements of the operators should not depend on parameter a chosen by this way. Let us mentally decrease 
the radius of a bound state, tq — > 0. Then the whole written scheme, as we have already noted, does not depend on 
a up to the values a = rQ. Therefore, the subspace H a can be identified with H due to the inequality R » a. In 
other words one can considers that "Pjj — > 1 for a — > 0. At the same time we do not break the quantum-mechanical 
description of bound states in virtue of the inequality a » rQ. From the physical point of view the inequality 
ro << R gives the stability domain for the bound states considered as elementary particles. The calculation of the 
following approximation is to be associated with accounting of the difference of the subspace H a from the space H. 
Finally we note that for an arbitrary vector | ) € H a the following evident relations are valid: 

C(x)|(x')| ) = 0, ( |C + (x)|+(x') = 0, |x-x'|<a, (12) 



C(x)|(x')5a - 0, i/ Q C + (x)e+(x') =0, |x- x'| < a, 
where £ and £ are any of the annihilation operators Xi,)(2, f]- 

III. THE STRUCTURE OF A(u, v) OPERATORS 

Here we consider the method for obtaining the operators A. Let the operator A represents a normal-ordered 
production of V»j(v), ^(u), (i — 1,2): 

i(u,v) =^+(u 1 )...^+(u 2 )^ 1 (v 1 )...^ 2 (v 2 )... (13) 

The operators of such type are the particle density operator /3j(x), 

Pi(x) = ?/>+(x)-0 4 (x), 

the momentum density operator 7Tj(x), 

i( l+( ^dU*) av3+(x) 

a Hamiltonian of the system, etc. 



The matrix element in the right-hand side of (|10f) may be written as the following vacuum average: 

(0|&(x x ) . . . ^ 2 (x 2 ) . . . <£ a (x) . . . i(u, v)tp+ (x) . . . V3 2 (x 2 ) . . . V3i(xi) . . . |0). 

Let us note that when calculating this average by using Wick's theorem, the quantity, which is averaged over the 
vacuum state has the meaning of a mixed T-product if to consider the operators ipi , . . . , ips, . . . , <p . . . referring to the 
moment of time +0, the operators ^ , . . . , ^ i ■ • • > ^ + • • ■ to —0, and the normal-ordered operator A(u, v) refereing 
to the moment of time 0. Thus, there is no need to place the contractions inside the expression for A(u, v). Let 

A>(y;y';u,v) = A b ( yi , . . . ,y 2 , . . . ,y, . . . ;yi, . . . ,y 2 , . . . ,y', . . . ;u,v) (14) 

be the analytic expression that corresponds to the diagram " 6" . For this diagram the operators with the arguments 
u are related to the operators with the arguments yi, . . . , y 2 , . . . , y, . . .. In addition, the latter arguments are spaced 
apart by the distances greater than a and coincide with some of the arguments Xi, . . . , x 2 , . . . , x, . . . ,. The similar 
statement should be also made concerning the arguments y' 1; . . . , y 2 , . . . , y', . • ■• Then, the operator A(u, v) acting in 
H is given, according to © , by 

I(u, v) = I RiR2RA b (y; y'; u, v^R^R', (15) 

b J 

where 

Ri = i[x + (yi)dyu R2 = Hx + (y2)dy 2 , R = l[fj + (y)dy, 

r'i = n*(yi) d yi' % = UWM> R' = Xlviy'W 

Here the summation is taken over all diagrams of the described type. 

If A(v) = 1 (see the proof of 0), then A(v) = 1 on the subspace H a . Let now A(u, v) = ^i(v). Then the only 
diagrams of the described type for the vacuum average (0| . . . ipi(y) • • • |0) are the following diagrams: 

A hx =^i(v)...^+(y / ), A b2 =^ 2 (y 2 )...^(v)...^+(y'). 

a a a b ba 

The expressions 

^ 1 (yi;v) = 5(v-y' 1 ), 

^6 2 (y2;y',v) = J dxidx 2 (p Q (xi - x 2 )S(y' - X)<5(v - xi)<5(x 2 - y 2 ) = (^(v, y 2 , y'), 

correspond to the above mentioned diagrams. Here <y9 Q (v, y 2 , y') is defined in accordance with (J2J. Therefore, 
according to l|15|) we have 



^i(v) = J 4yiAi(yi,v)xi(yi) + J dy 2 dy'A b2 (y 2 ;y',v)xt(y2)fl(y') = Xi(v) +Oi(v), (16) 

where 

Oi(v) = f dy<p(v,y)xi(y), £(xi,x 2 ) = <p a (x)f} a (X) (17) 



and X is given by J5J , x = xi — x 2 . Similarly, one finds 

4(v) = X2(v) + 2 (v), (18) 

where 

2 (v) = J dyxt(y)<P(y,v). 

When deriving l|16|). I|18|l . we have essentially used the inequalities (J5J. 

Now let us consider A(u, v) = ^j 1 " (u)V>i(v). In this case the following five diagrams 

Ai = $...$-(u)&(v)...$+ A 2 = Tp...i>...$f(u)ipi(v)...(p + ) A 3 = ^...ip+(u)Mv)...^ + ...<4> + , 

a a b b a b b c ca ab b c c a 

A 4 = tp . ..^(u)^i(v) A 5 =^...ip. ..^(u)i>i(v) . . .ip + . . .-0+ 

ab b c ca a bd d c ca b 

correspond to this operator (the indices 1 and 2 for -0 and ip + can be easy restored if to take into account © and the 
definition of contractions). The analytic expressions for these diagrams have the form 

Ai = Xi"(ui)xi(vi), A 2 = J dz 1 dz 2 <p(z 1 ,z 2 )5(z 1 - vi)xt (ui)%t ( z *) = J dz 2 £(vi, z 2 )xt (ui)xj(z 2 ), 
A 3 = J dzidz 2 ip + (zi,z 2 )(5(zi - Ui)x 2 (z2)xi(vi) = j dz 2 <£+(ui, z 2 )x 2 (z 2 )xi(vi), 
A-4 = J dzidz 2 dz\dz 2 (p + (zi, z 2 )ip(z[, z 2 )S(z 1 - ui)<5(vi - z' 1 )5(z 2 - z 2 ) = J dz 2 ip + (ui,z 2 )tp(vi,z 2 ), 
A 5 = - J dz x dz 2 dz\dz 2 (p + {z x , z 2 )ip(z\, z' 2 )5(z 1 - Ui)<5(vi - z[)xi (z' 2 )x 2 (z 2 ) 
= - J dz 2 dz 2 (^ + (ui, z 2 )<£(vi, z' 2 )xt (z' 2 )x 2 (z 2 ), 
whence, we find the operators (|15|l corresponding to the diagrams (|14p. 

k =Xi"(u 1 )xi(vi), i 2 =x^(u 1 )0 1 (v 1 ), I 3 = 0+(u 1 ) X i(v 1 ), i 4 + i 5 = di-(u 1 )(5 1 (v 1 ). 

In obtaining the latter expression, we have taken into account the anticommutative relations for x, x + ■ Next 
bearing in mind i|16[l. we obtain the final expression for A(u,v) = -0j t "(ui)^i(vi) that corresponds to 

i(u,v) =^+(u 1 )0 1 (v 1 ): 



Similarly, one gets 



^(ui)^i(vi) -» #-(ui)&(vi) = ^(ui)^i(vi). (19) 

-0l(Ul)-0l(vi) -> 0l(Ul)0l(vi) = ?/'l(Ul)'0l(vi), 

#-(tii)&(vi) -> 0+(u 1 )Vi 2 (v 1 )= ^+( Ul )4(vi), (20) 
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In general case the following formula is valid: 

K ) • • • $t (va ) ■ ■ • i>j m (v ira ) -> VA+ (ui, ) . . . 4 m (v Jm ) = it (ui) • • ■ <L (v im ) ■ (21) 

To explain this formula we note that each of i>j(vj) (or ^(u^)) entering A(u, v) (see l)13[l) is related to other operators 
■0 + (or ■0) not entering A(u, v) only by unique way, which leads to the binary relations (|19fl . 1|2U|) . Therefore, we come 
to l|21f) by sorting out all the operators ipi, tyv containing in A(u, v). 

The operators f/>j(x), ^(x') are anticommutative. For this reason there is a question concerning the consistency of 
(l2"0l . The anticommutativity of $i(x), $j(x') (and also Vv(x), %pf{x') ) represents the consistency condition of 

these formulas, 

{^(x)J,(x')} = {^+(x),v"+(x')} = 0. 

The validity of these anticommutative relations can be easily proved if to use the definitions (|16f) - H18|) f° r i>i ^ an d 
the commutative relations for x, x + and r) a , 17+. 

IV. THE OPERATORS OF BASIC PHYSICAL QUANTITIES 

In this section we consider the operators of basic physical quantities, which act in the Hilbert space H. Let us start 
from the density operator for particles of the first kind. The corresponding operator acting in the original Hilbert 
space H is of the form 

pi(x) = ?Aj f (x)'0i(x). 

Hence, in accordance with l)2(J|l . one finds 

~pi(x) = ^(x)^(x) = x^x)xi(x) + 6+(x)xr(x) + #(x)<3 x (x) + O+WO^x). 

Note that the operators with zero matrix elements in the subspace H a appear in the right-hand side of l|19|) because 
of points Ui and vi are close to each other. Since 

Oi"(x)xi(x)= / dy^ + (x,y)x 2 (y)xi(x), xlWOi(x) = / dyxf (x)x£(y)<p(x, y) 



and <p(x, y) differs from zero only for |x — y| < a, these operators according to l|12[l do not have the matrix elements 
in the subspace H a and therefore can be omitted. Using the permutation relation {x^(zi), X2(z 2 )} = S(z\ — z 2 ) the 
operator (5 1 ~(x)Oi(x) is written as 

6^ (x)Oi(x) = y dz<p + (x,z)fi(x, z) + y dzxdz 2 (p + (x, z 2 )(,c>(x, zi)xj(zi)x 2 (z 2 ). 

The matrix element of the second term is zero in the subspace H a because of Zi w z 2 ~ x (see dj). For this reason 
this term can be omitted. Thus, with the use of the method that was described in previous section we have 



pi(x) = xr(x)xi(x) + / dz^ + (x,z)<^(x,z). (22) 
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Similarly, if p 2 (x) = ^(xjifefx) represents the density operator for particles of the second kind, then /5 2 (x) — > /02(x), 
where 



P2(x) = x 2 f (x)x 2 (x) + J dz0 + (z,x)0(z,x). 



(23) 



Bearing in mind (|16|) and the assumption concerning the "small radius" of the bound state, one gets the following 
formulas: 



which allow us to obtain the densities operators for particles of the first and second kinds, 

Pi ( x ) = xt (x)xi (x) + »?+ (x)?) Q (x) , p 2 (x) = xt (x)X2 (x) + ?)+ (x)») Q (x) . 



(24) 



Thus, the operators 77+ (x), 7y Q (x) can be interpreted as the creation and annihilation operators of the bound states 
with quantum numbers a at the point x, and ^(x)?) a (x) as the density operator of the bound states. For example, 
the first formula from (|24|l have a simple physical meaning: the density of particles of the first kind is equal to the 
sum of densities for free particles of the same kind and bound states (each bound state contains one particle of the 
first kind). 

Consider a state vector $(X) = ?)+(X)|0), which specifies a compound particle at the point X (this state vector 
corresponds to continuous spectrum) . Then, in accordance with (|24(l , we have 

M 



($(X), / 5 1 (x 1 )$(X')) = <S(X-X') 



( M 
\m 2 



lf a I (Xi - X) 

1712 



For a wave packet 



the quantity 



Xn 



dX/ Xo (X)*(X), 



dX 



M 

(*Xo,Mxi)*Xo) = ( 

\m 2 



dX 



/x (X) 



/Xo(X) 



ip a I — (xi - X) 

mi 



should be treated as the probability density to find the first particle at the point xi if the atom is in a state \&x • If 
the bound state is localized near a point X (i.e., |/ Xo (X)| 2 — > <5(X — X )), then 



(*Xo,/5i(xi)*Xo) -» ( — 
m 2 



ip a I — (xi - X ) 

m,2 



Since (M/m2)(xi — Xo) = x = xi — x 2 , we come (as it should be) to the probability distribution for the space 
coordinate of the first particle in atom (the atom is at the point Xo). 

Let us find now the momentum density operator in the space H. In original Hilbert space, the momentum density 
operator 7Ti(x) for particles of the first kind is defined as 

* /„ n i ( l+i n^i(x) di>t(x) J 

-i(x) = --U (x)-^ ^^iW 
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Then according to (|20|l 



7Ti(x) -> 7Tl(x) = -- I Vl W— 7^ ^— Vl(x) 



Following the derivation of (|22H . (|23|) for pi(x) and /5 2 (x) we obtain 

ir l (x)=--f X i-(x)- s ax-Xi(x) -~ 2 dyU + (x,y)^ -—^y)], (25) 



~ I \ i f -+f \!Mx) ^X^(x) » , A i /" , dy(y,x) (9cp+(y,x) \ 

It is convenient for our further consideration to rewrite l|25|) . I|2tj|) in terms of the center of mass variables y = yi — y 2 
and Y = (miyi + TO 2 y 2 )/(r7ii + m 2 ): 



(x) 9xi"(x) „ \ i / ^vr/ v m 2 \ /-+/ vA 9<^(y,Y 



*i(x) = -dirW^ - ^T^Xi(x) - - dydYS I x — Y -fy J £ + (y,Y) 



2\ 9x Sx A v 7 2j ' " \ AT/ t w ' ' dy 



dy y M V 9Y 9Y 

^(x) = -- X 2 W— 53 5 X2(x) - - / dyrfYcS x - Y + — y <^ -<^+(y,Y)- 



dx 3x A ^ 7 2j J V M*J\ rv "'dy 



(27) 



dy r v ' Af V 9Y <9Y 

where </5(yi,y 2 ) = <£(y, Y). Note that in terms of these variables the operators pi(x) and p 2 (x) have the form 

fc(x) = x^(x)xi(x) + J dydYS (x - Y - ^y) 0+(y, Y)0(y, Y), 

fr(x) = X 2 + (x)x 2 (x) + / rfyrfY<5 (x - Y + ^y) 0+(y, Y)0(y, Y). (28) 
It is clear that (gEl 

can be expressed through the creation and annihilation operators ?)j(x), ?Jq(x) of atoms if 
to employ l|17fl . Taking into account l|27|l it is easy to find the operator n = 7Ti(x) + ^(x) of the total momentum 
density of the system in the approximation, where the radius of bound state is small, 

~ 1 f* + ( n 5 Xi(x) dxf(x) A \ A 9x 2 (x) dxt(x) - , A * /^+/ ^df) a (x) dfj+(x) „ 
7r = -2^ (X) ^ 9x- Xl(x) J " 2^ (X) ^X 9X-^ 2(X) J - 2^ (X) ^X 9x-^ W 

The third term in this formula is in accordance with interpretation of r)+ (x), fj a (x) as the creation and annihilation 

operators of the bound state with quantum numbers a at the point x. 

Finally let us consider a Hamiltonian in the space H . We suppose that this Hamiltonian has a standard form in 

the Hilbert space H and can be written as 

H = Ho + V, 

where H.q and V are the operators of kinetic energy and potential energy given by 

i 2 r 

V =2^ J rfxdx'^+(x)^+(x')^(x-x')^(x')^(x), 



<-j=i ' 
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and ^jj(x — x') is a potential energy of interaction for particles of the kinds i and j. After the similar calculations 
that led us to the expressions for p, 7r, we obtain 



2 



i=i 



2m,; \ / dx dx I ' ' <9x,; <9x. 



Next changing to the center of mass variables (see iJUJ) and noting that 

d d mi d d d m 2 d 

9x7 _ a^ + I79x' = ^dx" + aT ax' 



one gets 

2 



n o = ^^J rfx^- — — - - y dxdx^ (x)^(x)^- — — + — y dx— — — (si) 



where /i = m^/fmi + 7712) is a reduced mass. 

Let us find now V (V —*V). According to 1|2U[I we have 

2 



E / rfxdx'^ + (x)4 + (x')^(x-x')4(x')^(x) 



where V>;(x) = x«( x ) + Oi(x) (see (|2U|l '). Thus, can be represented in the form 

v = % + Vi + v 2 + v 3 + n, 

where Vk [k = 0, ...4) contains k multipliers of type x an d 4 — k multipliers of type O. The operators Oi(x) have, 
according to (|T%1) . the form 



Oi(x) = y dy<^(x,y)x+(y), 
01 (x, y) = <^(x, y), 2 (x, y) = <£(y, x), (32) 
where index i' is defined as 1' = 2, 2' = 1. Then Vq is of the form 
~ 1 2 f 

V = - 22 / dicidx 2 Uij(xi - x 2 ) 

x / dyidy2dy3dy4£+(xi,yi)(£+(x2,y 2 )£j(x2,y^^ (33) 



Note that the operators (p and in (|33|l are normally ordered, whereas x an( i X + are n °t normally ordered. Therefore, 
we put them in order using Wick's theorem: 

x*'(yi)xj'(y2)x^(y3)x^(y4) = : x*'(yi)xj'(y2)x/(y3)x 4 v(y4) : + : &'(yi)xi'(y2)x^(y3)x?(y4) = 

a a 

+ ■ x i '(yi)xi'(y2)x^(y 3 )x 4 v(y4) : + : Xi'(yi)xj'(y2)x^(y3)x?(y4) : + : x 2 '(yi)xj'(y2)xjt(y 3 )x?(y4) : 

a a a a a a 

+ ■ Xi'(yi)xj'(y2)xJ(y3)x^(y4) : + : Xi^yi)xj'(y2)xJ(y3)x-<(y4) : • (34) 

a b a b b {, a 
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The operator tp(pc,y) = (p a (x — y)??a( m 7 ^ : ^^ y ) differs from zero only for x sa y, (|x — y| < a). Thus, only those 
of ifi, for which |yi — y4| < a contribute to the integral over y in l|33|) . This means that the first term in (|34|) in 
virtue of (|12|l does not contribute to the matrix element of Vo taken between the states belonging to H a because 
of <£>i(xi, y4)xi(yi) < i ) = 0. Similarly, one can prove that the terms, which contain the single contractions in (|34|) do 
not give a contribution to the matrix element of Vo taken between the states in H a . The penultimate term in 1341) 
containing the double contractions does not also contribute to the above mentioned matrix element. Indeed, the 
penultimate term in (|34(l equals to <5(yi — y3)^(y 2 — Y4,)- In this case the nonzero matrix element exists for Xi w x 2 
and </3j(x2,y3)(/5i(xi,y4)<£> = in virtue of ljl2|l . Thus, only the latter term in (|34(l equal <5(y2 — y3)<5(yi — y4) can 
give the contribution to the matrix element of Vq. Therefore, not changing the matrix elements in H a , the operator 
Vo can be represented in the form 

V = - 22 / dxia'x2dyia'y2^ + (xi,yi)^+(x2,y2)^(xi - x 2 )<£j(x 2 ,y 2 )<£i(xi,yi), 
i,i=i J 

or, according to (|32|l . 

V ® = \J rfx i rfx 2rfyirfy2^ + (xi,yi)^ + (x 2 ,y 2 ) 

x |^n(xi - x 2 ) + z> 22 (yi - y 2 ) + ^i2(xi - y 2 ) + ^2i(yi - x 2 ) j<£(x 2 ,y 2 )<£(xi,yi). 

Similarly, noting that 

2 



Vj = \ ^ J dx 1 dx 2 dyvi j (xi - x 2 ) 

x \ ^(xi,y)Xi'(y)x J + (x2)Xj(x 2 )x l (xi) + ^(x 2 ,y)xi(xi)x^(y)xj(x 2 )xi(xi) + /i.c. \, 



^3 = r^ / dxidx 2 dyidy2dy 3 Ui j (xi - x 2 )| ^+(x x , yi)^+(x 2 , y 2 )^ i (x 2 , y 3 )x l '(yi)Xi'(y2)x^(y3)x»( x i) 

+ ^ + (xi,yi)^|(x 2 ,y2)^ i (xi,y 3 )x l '(yi)xi'(y2)Xi(x 2 )x 1 v(y3) + h.c. 



and performing the same derivation as for Vq, it is easy to verify (using the anticommutative relations for x, X + an d 
(13411 ) that we can consider V\ = V3 = not changing the matrix elements of V\ and V3 in the subspace H a . Next it 
is evident that the following formula is valid: 

1 ^ 

^4=9$^ / rfx i dx 2X l + (xi)Xj h (x 2 )^(xi -x 2 )xj(x 2 )xi(xi). (35) 

Finally let us find V2, 

V 2 = \ J2 I dx 1 dx 2 v ij {x l - x 2 )|6 t +(x 1 )6t(x 2 )x J -(x 2 )xi(x 1 ) + 6 l +(x 1 )x+(x 2 )6 J (x 2 )x i (x 1 ) + h.c. 
1,3=1 J 

+6+(x 1 )x+(x2)x,(x2)d i (x 1 ) + x+(xi)0+(x 2 )0,(x 2 )x 4 (xi) 
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It can be easily seen that the first two terms and the corresponding Hermitian conjugate terms do not give a contri- 
bution to the matrix element of V2 in the subspace H a . Therefore, we can consider 
2 



% = 2 ]C / dx i rfx 2rfyirfy2^(xi - x 2 )j ^(xi,y 1 )0i(xi,y 2 )xi'(yi)x/(x2)xj(x 2 )x-;(y2) 

+ t (x2 , yi)0j (x 2 , y2)x+ (xi)xjv (yi)xi (y2)xi (xi) 



The first and second terms in this expression contribute to the matrix element of V2 in iJ a under the following 
arrangements of contractions: 



: Xi'(yi)x|( x 2)Xi(x2)X?(y2) : + : Xi'(yi)x/(x2)Xj( x 2)x l t(y2) := ^v<5(yi-x 2 )^(y 2 -x 2 ) + (5(yi -y 2 )x J + (x 2 )Xj( x 2), 

a a b b a a 

: X l + (xi)Xi'(yi)X/(y2)x l (xi) := S( yi - y 2 )x+( x l)Xi( x l)- 

a a 

Thus, we have 

% = J dxidx 2 ^i 2 (xi - x 2 )(^ + (xi,x 2 )(^>(xi,x 2 ) + J dxidx 2 (^ + (x 2 ,y 2 )^(x 2 ,y 2 ) 

x |zyn(xi-x 2 )xi L (xi)xi(xi) + ^ 2 i(xi-y 2 )x5 f (xi)xi(xi) + ^ 22 (xi-y 2 )x 2 f (xi)x2(xi) + i'i 2 (xi-x 2 )x 2 f (xi)x2(xi)|. 

(36) 

The first term in this formula quadratic in field operators can be combined with the latter term in (|31() . As a result 
we obtain 

9v3*(x) dv? a (x) 



dxidx 2 i/i 2 (xi - x 2 )^> + (xi,x 2 )^(xi,x 2 ) - — J dxdX77+(X)7? Q ,(X) — |-— 



(9x 



| dxdXr)+(X)ry /3 (X)^(x)|-i-A x + z, 12 (x)}^(x). 



Since </?a(x) satisfies the Schroedinger equation 



-^A x + i/i 2 (x)|(^ /3 (x) = Epvpix), 
where ep are the atomic energy levels, the latter formula takes the form 

f dXY,e a fj+(X)f, a (X). 

J a 

Hence, taking into account 1)30(1 . 1(31(1 . ((36(1 . the Hamiltonian of the system TC is of the form 

/t — no + a int + /x int + n int , 

where 

j — 1 J a v 

is the Hamiltonian for free particles and bound states, and 
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T^lnt = J dx 1 dx 2 dy 2 'f + (x 2 ,y 2 )(p(x 2 ,y 2 )^ (vu(xi - x 2 ) + v 21 {x 1 - y 2 )Jx| (x 1 )xi(x 1 ) 

+ ( ^22(xi - y 2 ) + 1^12 (xi - x 2 ) )x 2 h (xi)x 2 (xi) (38) 



T<1nt = \ J dx 1 dx 2 dy 1 dy 2 ^ + (x 1 ,y 1 )^ + (x2,y2)^(x 2 ,y 2 )^(x 1 ,y 1 ) 

x <^ i/ u (xi - x 2 ) + j/ 22 (yi - y 2 ) + i^ia (xi - y 2 ) + 1*21 (yi - x 2 ) }, (39) 



^\nt = \ J dxidx 2 |z/ii(xi -x 2 )x^(xi)x^(x 2 )xi(x 2 )xi(xi) 

+ J/ 22 (xi - x 2 )x 2 f (xi)^ (x 2 )x 2 (x 2 )x 2 (xi) + 2z/ i2 (xi - x 2 )X] h (xi)xi(xi)x 2 h (x 2 )x 2 (x 2 ) i- (40) 



are the Hamiltonians of interaction. The Hamiltonian 7i l 1 „ t corresponds to scattering of particles of the first and 
second kinds by bound states; the Hamiltonian n\ nt corresponds to scattering of bound states by bound states; 
finally, the Hamiltonian Hf nt corresponds to scattering of particles of the first and second kinds by particles of the 
same kinds. The Hamiltonians of interaction (|38|l . I|39l) may be written through the creation 77+ (x) and annihilation 
f] a (x) operators of atoms by using of (|17|) . We want to emphasize that the obtained Hamiltonians of interaction do 
not lead to decay processes and formation of compound particles as it should be in the low-energy approximation. 
This fact reflects that atoms are absolutely stable in the main approximation. 

In conclusion of this section we address to the Galilean invariance of the developed theory. We suppose that the 
initial theory based only on the field operators \pi (x) , (x) , (i = 1,2) is Galilean invariant. This means that under 
the following transformations: 

&(x) -> #(x) = e^* vx ^(x) 

the operators for particles number Ni, momentum Vk and energy TL transform as follows: 

Ni — > N[ = Ni, i = l,2, 
Vk V' k =V k + v k (Ni +N 2 ), H^H' =H + v k V k + \m 1 v 2 N 1 + ^m 2 N 2 v 2 . 

Since i/'i( x ) an( i V'i( x ) meet the same commutation relations, they are related to each other by the unitary transfor- 
mation 

f j (x) = ^ i (x)D+ 1 

where the unitary operator U v , is defined by 8 : 



U v = cxp f — ivnii J dxxpi (x) — zvm 2 J dxx.p 2 (x 



(41) 
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Let us show that 



U v = cxp^— ivmi J ebcxp^x) — «vm 2 J dxx/5 2 (x)^ 



defines the Galilei transformations in the presence of bound states of particles (compare to l|410). To this end we note 
that this operator can be written, according to l|28(l . in the form 



U„ = 

where 



v = cxp^-iv^TOi J dxx+(x)xi(x) - ivs^j, 



s = mij dxidx 2 xi^ + (xi.x 2 )^(xi,X2) + m 2 J dxidx 2 xi<^ + (x 2 , Xi)<£(x 2 , Xi) 

= J dxidx 2 (mixi + m 2 x 2 )( y 5 + (xi,x 2 )(^(xi,x 2 ). 

Since 

<£(xi,x 2 ) = ^ tt (x)ry Q (X), y dxip* a (x)ip / a(x) = 5 a/3 , 
changing to the center of mass variables one finds, 

s=(to 1 +to 2 ) / dXX^^(X)i(X). 

Therefore, just as we expected 

Uy = expl -ivy^mj / dx£t(x)xt(x) - iv(mi + m 2 ) / dXX^ j)+(X)») a (X) J . 

^ i=l ^ ^ a ' 

The employment of this formula and the canonic permutation relations result in 

U v Ux)U? = e™ lVX X,(x) = x' t (x), U v r) a (X)Ut = e 4(mi+ " l2)vX ^(X) = j&(X). (42) 

Since the quantity <^ + (x, X)0(x, X) is not changed under such transformation, the Hamiltonians of interaction H} nt , 
Hi nt , T^int ( see CHt - (EDJ) are invariant with respect to (@2J. It is easy to verify that the Hamiltonian H for free 
particles and bound states transforms according to the law 

Hence, 

t/„W+ = H + v k V k + ~^-Nt + ~^-N 2 + — ——^ — — Nb, 

U v VkU+ = Vk + miv k Ni + m 2 VkN 2 + (mi + m 2 )v k N b , 

where Vk, N±, N 2 , N b are the momentum operator and the particle number operators (for free particles of the first 
and second kinds and their bound states). These formulas prove the Galilean invariance of the developed theory. 
Finally we note the validity of the following formula: 

U v Tri(x)U+ = 7Ti(x) + mivpi(x), 

which follows from ip7|). by using g2J|. 
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V. ELECTROMAGNETIC INTERACTION 

Here we consider the electromagnetic interaction assuming that the formation of bound states of particles is caused 
by Coulomb's (electromagnetic) forces. Therefore, the potential energy i/y(x— x') entering (|38[l - i|4U[) should be written 

as 

^.(x-x') = r ^^, i,j = 1,2. (43) 
|x — x'| 

We also introduce the additional interactions of particles (see below) with an external electromagnetic field A^ e ^(x, t), 
(p^ (x, t) and quantized electromagnetic field specified by a potential a(x) (Coulomb's gauge) 

A(x,i) =a(x)+A( £ )(x,i), 

where 8 



k A=l V 7 V 

(CkA is the annihilation operator of a photon with momentum k and polarization e^). As in section II, we identify 
the subspace H a to H assuming that the matrix elements of the operators in coordinate representation give the main 
contribution to the quantum electrodynamics processes and correspond to the space scale Ax > a (a — > 0, a >> ro). 
The density e(x) of a Hamiltonian in the second quantization representation has the well-known form 

e(x) = £/(x) + £ C (x) + ei#(x)^(x)<^ e ) (x, t) + ^ - ieM*, <)) ^ + (x) + ie.A^, *)) &(x), 

i=l i=l 1 ^ ' ^ ' 

where £/(x) is the energy density of a free electromagnetic field, £c( x ) is the energy density of Coulomb's interaction, 
the third term describes the interaction between the particles and external scalar potential. This formula leads to the 
following Hamiltonian of the system: 

H(t) =H Q + n int + V(t), H Q =H f +H p , (44) 

where 

H f = f dxe/fx) = ]>> fc (7+ A (7 kA , (45) 

^ 1 t 0#"(x)0&(x) 



i=l 



W int = / dxe c (x) (47) 

are the Hamiltonians for free photons, free particles, and Coulomb's interaction respectively. The operator V(t) 
represents a Hamiltonian that describes the interaction of particles with the electromagnetic fields A(x, i), ip^>(x, t) 

2 

V(t) = -J dxA(x,t)J(x,f) - i y dxA 2 (x,t)^ ^-6- !; (x) + y dx^ (e) (x,i)a, (48) 
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2 2 



J(x,t) = -A(x,t)V-*-fc(x)+j (x), jo(x)= x T-^7f l (x), (49) 

z — ' Hi; z — ' TO; 

i— 1 2—1 

where Ui{x) — e%pi{x) is the charge density for particles of the kind i, a = &i + <J2- This Hamiltonian is exact and 
describes (in this sense) the processes in which the bound states of particles are involved. Since a(x) commute among 
themselves and with V>i(x), \jjf{x), in order to find the effective Hamiltonian H(t) that describes electrodynamic 
processes at low-energies in the presence of bound states of particles we need, in accordance with section IV, to 
replace the operators a by "tpi, ijjf, a: 

H(t)=H(t)\ T , ., 

where 

^i(x) = Xi(x) +J dy^(x,y)xj(y), V> 2 (x) = Xa(x) + y dy^(y, x)£f (y), a = a 

(a(x) coincides with a(x) and acts in the space H). As a result the effective Hamiltonian is defined by (|44(l - 149|l . in 
which 

CTi(x) -> o-j(x) = e i/ 3i(x), 7Ti(x) -> 7Ti(x), ^i(x) -> ^i(x), (50) 

moreover, the operators /3j(x), 7Tj(x), V'i(x) containing the creation 77+ and annihilation 7/ a operators of bound states 
are determined bv l|27 |l . QR fy . (P|) - (Pjl . 

Let us define the operator of magnetic field 

H(x,t)=h(x)+H( e '(x,t), (51) 

where 

h(x) = rot a(x) , H (e) (x, t) = rot A (e) (x, t) , 

and the operator of electric field 

E(x,t) = e(x)+E^(x,t), 

where 

e(x) = e f (x) +e'(x), EW(x,t) = -^-A^(x,t) - 

at ox 

The transverse e*(x) and longitudinal e'(x) components of the electric field created by particles are defined as follows: 



e'(x) = -^Uo(x), e 4 (x) = -a(x) = -*[£(<), a(x)]. (52) 
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Here ao(x) represents the operator of the scalar potential, 



ao(x) = j /rfx R — -. 

Next it is easy to verify that 

[r) a (X), 77+ (X')] - 5 a0 S(X - X'), ^W^( x ') = ^( x ~ x ') 

a 

lead to the following commutation relation: 

[0(x, X), 0+(x', X')] = 6(x - x')<5(X - X), 



(53) 



(54) 



where (p(x, X) = <^(xi,X2) is determined by l|17|) . Hereinafter we assume that a set of wave functions ip a (x') is 
complete. Therefore this set accounts for the wave functions of continuous spectrum. In addition, if we introduce the 
operator 

n(x,X) = £ + (x,X)<£(x,X), 

then 

[n(x,X),n(x',X')] = 0. 

Thus, it is clear that the Hamiltonian Hint, which describes the Coulomb interaction between particles and bound 
states meets the commutation relation 

[Wi„{,/3j(x)] = 0. 

Let us find now a commutator [p.j(x),7io] (* = 1,2), where f>i(x) and 7t!o are defined by Ij28(l . 145(1 . 1(46(1 . In order to 
do this we consider p\{x) in the form 

3i(x) = ^ 1 (x)+3?(x), 

where 

^(x) = X + (x)xi(x), ft'(x) = |dydY ( 5(x-Y-^y)^+(y,Y)^(y,Y). 

Then in accordance with the commutation relations for the field operators Xi(x), 

fc(x),X+(x')] = ^(x-x'), 



one gets 



^(x), / dx 



2 



2m 7 - dx! dx' 

3=1 J 



2m\ dxk 



dxk dxk 
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Next according to (|30|l we have 



«x), Wo] = - y ^ (x - Y - ^yj — ^ + (y, Y)^^ + -^^(y, Y) 

+ 2M 7 dydY5 ( X - Y - M y J 5n ^ (y ' Y) —dYT + dY k 
Bearing in mind 

dy k V M 3 ) M dx k \ M J dY k V M J dx k \ M. 

and taking into account the definitions of pi, iti (i = 1, 2), one finds 

i[pi(x),W +Hint] = — divfri(x). (55) 

771l 



Similarly, 



z[/3 2 (x),W + Hint] = — div7r 2 (x). (56) 
ra 2 



Let us find now a commutator [pi(x'), -fr k (x)], (i, k = 1, 2). With the use of (|2Tj) we have 



Noting also that 



we get finally 



[^(x'),^(x)]=^(x)^5(x-x'). 



^(x'),^(x)]=^(x)|-5(x-x'), 



5 

[p 4 (x'),-7r fc (x)] = iS ik pi(x)— £(x- X ). 

ax 



These commutation relations along with l|49|l result in 

fio(x),^(x')] = -iJ2—*i(x')^-6( X - x'). 

' ?77i OX 

2—1 

Thus, according to l(56^l we have 

£(x) = i[W(t),£(x)] = divJ(x,t). 
The definitions 1)52)) . 153[) of the longitudinal component of the electric field along with the above formula give 



e'(x) = i[W(t),e'(x)] = — div y ^'j^rr- (57) 

Let us find the equation of motion for e*. The use of commutation relations for creation and annihilation C^a 
operators of photons, 

[CkA, C'k'A'] — ^kk'^AA' 
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lead to 



e*(x) = -i[#(f),&(x)] =*EE(^T i ) (e^C'kAe^-ft.c.Y (58) 



k A=l 

4*(x) = -t[W(t), e*(x)] = rot h(x) + i[V(t), e*(x)], 

and also 

[a 4 (x),e*(x')] = -^^x-xO-ift^i^j. 
This formula and the definition of V(t) allows us to find the commutator: 

i[V (t),6*(x)] = -47rl(x,t) - |-div / dx'^4. 

ax 7 l x — x I 

Therefore, it follows from l|57)l. (JHEJl that 

e(x) = roth(x) -4ttJ(x,<), h(x) = i[W(f),h(x)] = i[H f ,h(x)}. (59) 

Upon calculating this commutator with the use of (|51|l . one obtains 

h(x) = -rote(x). (60) 

In addition, f5l|l . fB^ give divh(x) = 0, dive(x) = 47rcr(x). These equations along with represent the 

Maxwell-Lorentz equations for the operators of electromagnetic field 

2e(x) = roth(x) - 47rJ(x, i), dive(x) = 47ro-(x), 
h(x) = — rote(x), divh(x) = 0. 

Notice that gSJ, lEU>, © result in 



2 

k A=l J ^ 



where we neglect the unessential term J^kA 

Up to now we used the Schroedinger representation, in which the Schroedinger equation has the form 

In this representation a dot above the certain operator fc(x) (see ||S5J) does not mean the differentiation with respect 
to time. It means the following: 

S(x)=i[%(x') J t),S(x)]. (61) 

Here the Hamiltonian TL depends on dynamic variables x(x') (these variables unite x an d v) an d time t. Introducing 
the unitary operator S(t,0), 

= #(*(*'), o), 5(o,o) = i, 
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we can write the solution of the equation for <I>(i) in the form, 

$(f) = S(t,0)$(0). 

In the Heisenberg representation the vector $(0) = iff is taken as a time-independent state vector. Then, the time 
dependent Heisenberg operators 

satisfy the equation 

i^ = [x(x),fi(x(x),t)}. (62) 
Since the external fields H^ e \ E^ e ' satisfy the Maxwell equations, 



at 



= -rotE (e) , divH (e) =0, 



= rot H (e) - 47rJ (e) , div E (e) = 47R7 (e) , 



at 

the operators of total fields E, H in the Heisenberg representation satisfy, in accordance with (JfjJJ), the equations 

ati 

— = -rotE, divH = 0, 

at 

— = rotH-47r(j + J (e) ), divE = 4tt(& + cr (e) ), (63) 

where J*- 6 ) and are the extrinsic current and charge densities. 

Let us write out the expressions for the charge density operator <x and current density operator J in terms of 
creation and annihilation operators of particles and bound states. Bearing in mind 1|50|) we have 

2 



0"( X ) = X^'W' 

i=l 

2 2 2 

5(x) = -A(x) — Mx) + £ — **(x), (64) 

' — ' TUi * — ' TO,; 



where the particle density operators Pi(x) are defined by l|28|l and the momentum density operators 7Tj(x) by J27J. 
The first terms in expressions for /3j(x) and 7Tj(x) define the contribution of particles, whereas the second terms define 
the contribution of bound states (see l|27|) . (|28|l ). 

Let us find now the equations of motion for the field operators of particles Xi( x ) an d bound states 0(x,y). In 
the Heisenberg representation these equations are obtained from the general formulas l|62|) . Indeed, it follows from 
©-gDJl that 

Ki(x),W; n t(x(x'))] = |y dx 2 ^n(x - x 2 )/5i(x 2 ) + y dx 2 ^i 2 (x - x 2 ) j 02(x2)|xi( x ) 5 
[X2(x),7i mt (x(x'))] = |y dxii/ 22 (x- xi)/5 2 (xi) + J dx 2 i/ 2 i(x - Xi)pi(xi) jx 2 (x), 

[<£(x, y), Hint (X(x'))] = ^ rfX2 ^n( X _ X 2) + ^12(y - x 2^ Al(x 2 ) + 



y dx 2 ^ 22 (y - x 2 ) + i/ 2 i (x - x 2)^/02( x 2)|<^(x,y). 
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These formulas allow to find the equations of motion I162II for Xi( x ); X2(x) in the Heisenberg representation, 
«ii(x,t) = j-2~(J~: -*eiA(x,t)^ +ei^ (e) (x,i) +£>i(x,t)|xi(x,<), 



{~2^ ~ *e 2 A(x, i)^ + e 2 ^ (e )(x, t) + 2 (x, t) }x2(x, i). 



It is easy to find also the equation of motion for the field operators </?(x, y) of bound states, 

i<£(xi,x 2 ,*) = ^£ Q ?7 Q (X,i)<p Q (x) + {-^— ( ^ ieiA(xi,t) J - — !— ( ie 2 A(x 2 ,i) 



2 



(65) 



2m! \9x! ' J 2to 2 V<9x 2 

+ ei<p (e) (xx , t) + e 2 <p {e) (x 2 , t) + u(xi, x 2 , t) |^(xi , x 2 , t). 

Here 

Oi(x,i) = y dx 2 iALi(x- x 2 )/3i(x 2 ,i) + y dx 2 i/ 12 (x- x 2 )p 2 (x 2 ,i), 
w 2 (x,t) = y dx 1 i/ 22 (x-x 1 ) / 3 2 (x 1 ,i) + y dx X i/ 21 (x- x 1 )p 1 (x 1 ,i), 
u(xi,x 2 ,t) = J dx'^v n (xi - x') + ^i 2 (x 2 - x')^pi(x',t) + y dx'^ 22 (x 2 - x') + 1*21 (xi - x')j/S 2 (x',£). 

The operators ?) a (X,t) related to <£(xi,x 2 ,t) by 0(xi,x 2 ,t) = ( < 9 ) g(x)fy ) g(X, i). 

The equations of motion 1631) . (|65l) are gauge-invariant. Indeed, under the following transformations of fields Xi, 

Xi(x,t) - xi(x,t) - e ieia ^^Xi(x,t), fc(x,t) - ^ 2 (x,i) = e 4e2tt < x ^x 2 (x, i), 

0(x, y, t) -> £'(x, y, t) - e <eia(*.*)+i«a(y,t)^( X) yj ^ (66) 



and electromagnetic field, 



A(x, t) -> A'(x, t) = A(x, t) - |-a(x, i), ^ (e) (x, i) - ^ (e) '(x, i) = ^ (e) (x) - 

ox at 

the operators /3j(x) are invariant, 

&(x,t) -> /3-(x,f) = pi{x,t), 
whereas 7Ti(x) transform (see ||55J)), by the law 

. a 

7Tj(x,t) -> 7Tj(x,t) = 7Tj(x,t) + e;pi(x, i) — a(x, t) . 

ax 

As a result we have (sec (jHU) 

J(x,i) -»■ J'(x,i) = J(x,t), H^H' =H 
that provides a gauge invariance of the Maxwell-Lorentz equations l|63|l . Next noting that 

A - le! A(x, t)\ fc(x, t) - f |- - iei A'(x, t)) #(x, t) - e^ x <*> f A - 8 e t A(x, fc(x, *), 



ax / \ ax / \ ax 
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we make sure of gauge- invariance of the equations of motion for the field operators of particles Xi ( x ) • Taking into 
account the transformation law for the field operators of bound states <p(x,y), it is easy to prove a gauge invariance 
of the equation of motion for <£(x,y). 

Since tp(xx,x.2,t) = tp a (x)fj a (X.,t), we have under the gauge transformation l|66|) . 

where xi = X + (777,2 /M)x, X2 = X — (mi/M)x. Hence, we come to the transformation law for the operators r/ a (X): 

f)a(X,t) - f,' a (X,t) =^L afl (X,t)^(X,t), L Q/3 (X,i) = / dxv *( x ) e <«iapt+("«/M)x,t]4« a Bpt-( mi / J lOK 1 t] ¥ , /j(x)i 

/3 J 

In conclusion of this section let us note that "0i( x ) transform under the gauge transformation as well as x«( x ) : 

^(x,i)-^(x,t)=e le ' Q(x ' t) ^(x,t). 



VI. THE PROCESS OF SPONTANEOUS RADIATION BY ATOM 



In this section we employ the obtained formulas in order to find the probability distribution of spontaneous radiation 
of an excited atom. To this end we use the general formulas for the scattering matrix S in terms of the T-operator: 



/oo 
dE8{E - H )f {+) {E)5{E - H Q ), 
-OO 



f {+) (E)= lim f(E + irf), (67) 
77— >+o 



where T(z) is found from the equation 



f(z) = V + VRof(z), R = — (68) 

z - Ho 

In the first order of the perturbative approach over charge we have 

(/|f(*)|i> « (f\V\i), 

where 

K)=^ P |o), \f) = c+ x fj+ p \o). 

Here | i) is the initial state (in which atom is in the state a, p) and |/) is the final state (in which atom is in the state 
a',p' and photon is in the state A, k). According to l|48|) . I|5U|I we can take V in the following form: 



V = -j dxa(x)J(x), 

where the current of bound states is defined by 



2G 



By using of |@BJ , (EOJ w e find 

1/2 



The formulas JEJ, lead to 

j(k)= / dxe- 4kx j(x) = ^ X! ^ 1P1 ^2P2 fga lQ2 (k) + Pl 2 ^ P2 gctia2 (k) j A( Pl - p 2 + k) 

rti rti rvT rvi ^ ' 



Pi ,P2 "1 ,C*2 



and, consequently, 

1/2 



(/l^> = - (^) 4 A) * (g aa ,(k) + ^j^wflo) A(p - p' - k), k = p P ', 
where A(p) = <5 p0 , and 

g a p(k)= J dy(^e 1 e ik ^ m2 / M )+e2e- lk ^ mi / M ^;(y) W (y), 

** 00 = 4 / dy(^e^(-/^ - ^e"^ W)) (<(y)^ - ^W)) ■ (70) 

The matrix element of S- matrix is defined, in accordance with (|67|) . by 

(f\S\i) = -2iri{f\P\i)5(e(p) - e(p') - 

Next using the standard method we obtain the following expression for the probability of transition per unit time 
from the initial state ?7q P |0) into the final state flp p '^kx\^) : 



S(uj k + £/3(p') - £ a (p))dk. 



In the domain of small k (kro << 1, see section II) the Fourier-components of g Q( 3(x), g a p(pi) can be written in the 
form 

g«(i(0) =(—-—)[ dy<P* a (y)&Pfi(y) =(—-—) <a|p|/3>, ff«/j(0) = (ex + e 2 )8 aP . (71) 
\mi m 2 J J \mi m 2 / 

If e 2 = —e\ = — e (this case corresponds to the hydrogen atom), then <7 Q/ 3(0) = 0, g Q ^j(0) = (e//i)(a|p|/3) . The 
Schroedinger equation for the wave function ip a (y) gives 

HpI/3) = w(£a - ep) J dy(p* a (y)y(p p (y). 

Therefore, 

ga/?(0) = i(e a - ep)(a\d\P), 
where d = ey is the dipole moment of the atom. Hence, 

dw^ f = 27ru k \(a\e£ ) d\f3)\ 2 5(Lu k +e p (p') -e Q (p))dk. 
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Since e a (p) — e a +p 2 /2M, we have for the atom of infinite mass (7712 — > 00, M — > oo) 5 : 

dwi-^f = — \e^g a0 (k)\ 2 5(uj k +ep- e a )dk, 

where 

g ^ (k) = "2^T J ye V-M-dj ar^M 

(mi is the mass of light fermion). 

VII. SCATTERING OF PHOTONS AND FERMIONS BY ATOMS 

Here we study the scattering process of long- wave photons by atoms. For this purpose let us use the general formulas 
(lfT7Jl . (JSHJ, where the Hamiltonian of free particles is defined by 

2 2 

#o=£ / rfxVx+(x)Vx 4 (x)+^^ Wfc C+(7 kA + ^ / d J c(^-V^+(x)V^(x)+e a ^+(x)^ a (x)Y 

j=l J k A=l a J \ m ' 

and the Hamiltonian associated with electromagnetic interaction by 

V = V' + V". 

In accordance with l|48[) . (J5DJ, 

V' = \( dxa 2 (x) f — 3 x (x) + — ^(x)) , f " = - / dxa(x)J(x), 
2 7 \ m i m 2 J 7 

j(x) = — *i(x) + — * 2 (x) (72) 

(further we assume e\ = — e 2 = e). Since V' ~ e 2 , V 7 ' ~ e, in the second order in charge e, the T-matrix can be 
written as 

t{z) = V' + i?"Mz)V", r q = 1 



z - Ho 

(within the considered approximation it is not necessary to take into account the fields Xi( x )). Bearing in mind J75J, 
()28p. the operator V"' can be written as 



Pl,P2 

where 



«a/?(k)= / dy{— e lky(m2/M) + — e- lk ^ mi / M )j^(y) W (y). 

Similarly, 

t>" = -^ £ ^(Pi)^(P2)W(Pi+P2,Pi-P2) J dxaWe"^-^, 



Pl,P2 
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where 

W(Pi +p 2 ,Pi - Pa) = g a p{pi - P2) + Pl 2 ^ f P2 .9a/3(Pi " Pa). Pi-P2=k, pi + p 2 = P, 

moreover, g Q( g(k), g a p(lt) are defined by (j7U|) . 

Let C^xVa (p)|0) and (?k<A'^a'(p')|0) be the initial and final states of the system, 

\i)=CZ x f)i(p)\0), \i) = C£, x ,fj+(p>)\0). 

Then the matrix element (f\T(z)\i) can be represented in the form 

(f\f(z)\i) = {f\H) + {fft''Mz)t ,, \i), 

where 

(f\H) = —-^^^(p + k - p' - k'), 

~ - ~ 47T 2 ?' 1 

(/|V"i2o(z)V"|i) = -— - -=i& a A(p + k p' k'), A(k) = 5 k0 , 
and the matrices R',„, R'' ln , are given by 

<-q = 4 A)e k A) W(k), k = p' - p, 

ff „ _ v / (el^(p / ,kQ)*(eI^(p,k)) (e'I a/3 (p,-k'))*(eI et , ;3 (p / ,-k)) | 
Q ' a y\ io + e a (p)-s p (p + k) + _^ + £a (p)_ £/3 ( p _ k /) J- 

Here we have taken into account that u> + e a (p) = 10' + £p(p') is a consequence of (|67|l and assumed E = u> + e a (p). 
According to our approximation fcro << f (A >> ro) we have 

e 2 

qapQt) W — £ a/ 3, W 0, 



Mk)---/^(y)^ M = f(«lPl 



whence 

2 

*&a« -e^e^''*^, I^(P,k)~g Q/3 (k)«-(a|p|/?), (73) 

where p = — iV is the momentum operator. Noting that e a (p) = e a + (p 2 /2M) we neglect the kinetic energy in 
denominators of R n a i a (the kinetic energy is small in comparison to binding energy of bound states; see section II). 
The result is 

R " = e " V ( (<x'\e'*p\P)(0\ep\a) , (a'\ep\l3)(l3\e'*p\a) \ 
a ' a M 2 p I u + e a -ep -uj> + e a -e f3 J ' 

Let us express this formula through the matrix elements of the dipole moment of atom. It is easy to verify (by using 
of dJ that 



"e'*p|/3)(/3|ep|a) 2 / _ . ww' 



(a'|ep|/?)(/3|e'*p|a 



-a/ + e a - ep 



f i' 2 [e a -ep-L J ' + ) (a'|e'*y|/3)(/3|ey|a), 

fJ 2 (e a < -e p +u/ + ) ( a '|ey|/3)</3|e'*y|a) 

V -u' + £ a -EpJ 
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Therefore, the quantity R'^i a is of the form 



(a'|ed|/J)(/?|e'*d|a) (a'|e'*d|/3) (/3|ed| 



-lu' +e a - Bp 



where d = ey is the dipole moment of atom, and 

Qa 

Next, using that 



E|( e « - £ /3 - to')(a'\e'*y\0){P\ey\a) + + e a , - e /J )(a / |ey|/3)</3|e , *y|a)|. 



ye, 



2// 



-ep, 



one finds 



According to (|73|1 . the quantity -R^'a is equal to R' a , a — ^e'*eS aa i . Thus, the matrix R = R' + R" takes the form 

' (a'|ed|/?)(/3|e'*d|a) (a'\e'*d\(3)((3\ed\a) 1 



Ra'a = WUl' 2^ 



-U)' + S a - Ep 



u) + e a - ep J 

and the amplitude of transition from the initial state into the final state is given by (see 5 ) 

47T 2 i 1 



Si 



Ra' a S(e a (p) + u}-e a ,(p')-uj')A(p + \s.-p' -k') 



V \/ljJU) 

The differential scattering cross-section is of the form 

dai-,f = — !— \R a ' a \ 2 S (e a (p) +LJ- e Q '(p') - w') dk'. 

Let us consider the scattering of an elementary particle by compound particle (the scattering of an electron by 
atom) . Due to the structure of the Hamiltonian of interaction l|38l) it is sufficient to study the scattering of particle of 
the first kind by bound state. The matrix element of S-matrix in the first nonvanishing approximation over interaction 
is defined by 

S i ^ f = -2m{f\k\ nt \i)6{e i -e f ), 
where H* n4 is given by l|38f) and consequently, 

(f\k nt \i) = ^J dx 1 dx 2 dy V , Q (x 2 - y)e l P x ^,(x 2 - y) e -*P' x L n (x x - x 2 ) - ^(xj - y^e**^-*). 
Here (/| = (0|?7 Q '(p / )xi(k' 1 ) and \i) = xf (ki)fj+ (p)\0) . The simple calculations lead to 

(f\AlnS) = ^A(p + k- p'- k') J dx^*,(x)^(x) {ni(q)e- iqx(W2/M) + ^ 21 (q) e ' iqx(mi/M) } , q = k' k. 
For Coulomb's interaction of particles we have 



"ii (q) = -"2i (q) = "(q) 



Aire 2 
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Therefore, in the long- wave approximation (qro << 1), one gets 



(f\Hint\i) = 4^e^A(p + k - p' - k')^qd a , a , (74) 



where d = ex is the dipole moment of atom (d a > a = J dx< / 9*,x(x)( / 9 ct (x)). Note that the Hamiltonian 

V = J dxidx 2 x + (xi)x(xi)iV Q (xi - x 2 )7}+,(x 2 )7) Q (x 2 ), 

which describes the interaction between particles and bound states leads to the same result. Here x, f\ a are the 
field operators of particles and bound states respectively, and v a i a (x) — (e/x 3 )xd Q ,/ Q . The relation (|74|l defines the 
interaction energy of the dipole moment d and electric charge e, which are at the distance x from each other. 

VIII. THE VAN DER WAALS FORCES 

In this section we investigate the forces acting between neutral atoms being in the ground state (the van der Waals 
forces) on the basis of the developed formalism. In order to solve this problem let us address to the Schroedinger 
equation that determines the energy spectrum of the system, 

7i$ = E$>, H = H a + V. (75) 

Here Ho and V — H\ nt + H\ nt + Hf nt are defined by lj3*7)) . (|3"8|) - (|4*U|> . Since the system in question consists of two 
atoms, we should seek the solution of Q75JI in the form 

$ Q/3 (X,X') = £ / dYdY'K a0;Xp (X,X'; Y, Y')r}+(Y)r)+(Y')<f> . (76) 

The Hamiltonian of interaction V is equal, in accordance with l|38|l - (|40|) . to 

V = \ J dX^77+(X)j7+(Y)7) 7 (Y)%(X)G, 7;a ^(X - Y), (77) 



where 



Ga 7;ct /3(X- Y) = J dxdy^*(x)^(y)^ 7 (y)^(x)|i/i2 f X - Y - 



TOix + m 2 y\ /„ , miy + m 2 x 
)+u 21 X-Y + 



+ vn (x-Y + ^(x-y)) +i , 22 (x-Y-^(x-y 

We suppose that the kinetic energy of atoms is small in comparison to the energy levels |£ Q | (e a < 0). In this case, 
according to l|37fl . the operator Ho can be represented in the form 

Wo = ]T / dXe a fi+(X)fiJX). 

a J 

It can be easily seen that 

4ry+(Z)77 p +(Z')$o = (e x + £p )77+(Z)ry p +(Z')$ , (78) 
ff;+(Z) ? ) p +(Z')$o=E G V; Q /3(Z-Z')^(Z)^(Z')<i>o. 

a/3 
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These formulas show that we can seek the solution of <|75ll in a more simple, rather than l|76|l . form 

<J> Q/3 (X, X') = £ K aP , Xp (X, X')»£( x W x/ )*o, 

\ P 

so that the coordinates of atoms have the definite values in the state $ Q ^(X, X'). Upon substituting this expression 
into (|75|l and using i|78|). one finds 

K a p- n s{£~t + £ s) + '^ j K a p.\ p G\ p - 1 s(Z — 71) = E a pK a p.^s- 

Xp 

A perturbative approach for this equation can be easily developed in the domain of great | Z — Z'\ when the quantity 
G\ pn s{Z — Z') becomes small (see lj77)l). Expanding -ftT Q /3 ;7 5 in G, 

K a p- n s = Kl tl3 . 1& + K a p iyg + K a/3 . lS + ...., 
E a p = -E^ + -E a/3 + E a(3 + .. ., 

Ka/3;-fs( £ 1 + £ s) = E®pK°p. lS , 



one obtains in the zeroth order 



whence 



Ka/3;y8 ~ K°pS aJ Sps, E®p — E a + Bp. 

Taking into account this result we have in the first approximation, 

K a^ns( £ i + £ s) + KapGafrysiZ - Z') = (e a + sp)Klp. lS + E^pK^pS^Sps- 
By setting here a = 7, /? = 8, one gets 

Eip = G a p., a p(Z - Z'), (79) 

and for a, (3 7^ 7, (5, 

i^ = i^ , ^(Z-ZO 0,^7,* (80) 

"TO^ Bry tS ) 

The second order of the perturbative approach give 



Kaposi 6 ! + £ s) + X] K af};XpGxp^8(Z - Z') - (e Q + Ep)K%,p. n s + E aP K a^5 + E lpK°p n 
Xp 



Taking here a = 7, /3 = 5, one finds 



K a(3;XpGxp;af3{'Z' ~ Z') - E X a pK a p. a p + K^pE^p, 

Xp 



whence, according to J75J), i|5U)) . we have 



P 2 _ Y^'G a /3;Ap(Z - Z')G\ p;a p(Z - Z') 

^ " £ ( £Q +e,- £ A- £p ) ' (81) 
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The prime above the sum means that the terms with A = a, p — [3 are omitted. The state vector $ Q ^(X, X') in the 
main approximation of the perturbative approach is determined by 



$ a/3 (X,X') = ^+(X)7)+(X')*. 



(the constant K®p may be found from the normalization relation, ($> a p,$> a /3) = !)■ Formulas (|81|) for E^* and 
E?g give us the corrections to the energy levels E^g — e a + ep. It follows from the obtained formulas that the energy 
of two atoms being in the ground state a and spaced apart for sufficiently long distances is defined by 

E aa — 2£ Q + G aa - aa (Z — Z') + ^ — aa,\p{ ) — Xp,aa( _ _|_ _ _ _ (82) 

Ap 

Let us prove now that G aa - aa (Z — Z') = 0. In doing so we use the following formula 9 : 

oc / \ n 



\/i? 2 - 2i?/ra + p 2 |R-p| i? 



where x = cos$, $ is the angle between vectors R and p, and V n (x) are the Legendre polynomials. Noting that (see 

G3) 



G QQ;QQ (Z - Z') = y dxdy|^ Q (x)| 2 |^(y)| 2 |^ 12 (z - z' - 



mix + m 2 y\ / . miy + m 2 x 

1 ^21 Z - z 



M / V M 



+ Vn (z-z' + ^(x-y))+. 22 (Z-z'-^(x-y))}, 
we come, taking into account the spherical symmetry of |<^ a (x)| 2 , to 

G aa;aa (Z — Z ) = 0. (83) 

We have also employed here that according to l)43|) 

u ab (Z ~Z) = , z _ z ,, , ei = -e 2 = e. 

Let us find now G QQ;Ap (Z — Z'), 

G QQ;Ap (Z - Z') = J dxrfy^(x)^( y )^(y)^(x)|^ 12 (z - z' - m ^ x + m ^ + v%x (z - z' + miy + W2X 



, n (Z-z' + ^(x-y))+. 22 (z-z'-^(x-y))}. 



The presence of multipliers (p a {y) and </3 Q (x) makes it possible to expand the expression in braces in powers of x, y. 
As a result we have 

G Q q ; /3a(Z - Z') = ^ 3^7j3 ^-3(nd ( g Q )(nd AQ ) + (d^dAa) 

where 



Z — Z' f 
n = \ z _ z ,y d fja =e J dxx.<p*p(x)tp a (x). 
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(d^ Q are the matrix elements of the dipole moment of atom). Thus, the potential energy of interaction between atoms, 
V(Z - Z') is defined by 6 

V(Z-Z) = E aa - 2e a = jy—yTwl^ " 



|Z-Z'|6^ 2e a -e p -e x 



< 0. 



Since ep,e\ > e a , V(Z — Z') < 0. Therefore, the attractive forces (the van der Waals forces) act between the neutral 
atoms at long distances. 

Consider now the scattering of a compound particle by compound particle, which are in the ground state. We 
start from the scattering matrix representation in terms of T-operators l|t)7|) . Accurate to the terms, which are of the 
second order in interaction V between the atoms (see (jHEJ), one gets 

/oo ^ ^ />oo ~ -i ^ ^ 

dES(E -H )VS(E -H ) - 2ni / dES(E-H )V = V5(E-H ). (84) 
-oo "'-oo E — Ho + ii] 

Neglecting the kinetic energy of atoms we replace in the resolvent (E — Ho + ii])^ 1 the Hamiltonian of free particles 
Ho by the Hamiltonian Ho of atoms in the rest (see H37(l 'l 



(in the low-energy limit we can not do this under the argument of (5-function). 

Since we study the collision of two atoms and the number of atoms is conserved during the process of collision, 
we can restrict ourselves by considering of this process in the two-particle subspace with the following completeness 
condition: 

i J d Xl dx 2 J2 ^ 1 (^i)^(x 2 )|0)(0|r) Q2 (x 2 )r) Ql (x 1 ) = 1. (85) 

Ql,Q2 

In this case formulas l|78|l are true. Since the atoms are in the ground state, l|83|) is also valid and, therefore, the 
matrix element of S'-matrix turns to zero in the first approximation of the perturbation theory. 

The second approximation of the perturbation theory is defined, according to l|84|) . I|85|) . by the formula 

V F l,J = \Y, £ £ / i~ - s a2 + i, Ga '^ft (xi - X2) 

E-Ho + irj ai,a 2 ft.ftj 71,72 J 1 2 ' 

x ^/3 1 ( x i)^/3 2 ( x 2)|0)(0|r) 7l (x 1 )r} 72 (x 2 )G+ 72 . QlQ2 (xi -x 2 ). 
In the two-particle subspace this formula is equivalent to 

V i V= ( dxidx 2 7 ?/3i( x i)^/3 2 ( x 2) t/ ft/52;7i72( x i - x 2;J E;)^ 7l ( X i)?7 72 (x 2 ), 

E-Ho+if] J ft,^ 71,72 

where 

^1&;7172( X 1 -*2]E) = Y tti — T~ G 7i72:aia2( X l _ X 2 )G Qia2;i g l/ 3 2 (x X -X 2 ). 

ai a 2 ai a2 + ' 

Since the initial and final states of a particle are the ground states, the scattering matrix in the two-particle subspace 
can be represented in the form 

r oo 



/oo 
dE6(E - H )VeffS(E - Ho), 
-OO 
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where 

Vef/ = ~ J dXi J rfX 2 f?+ (x 1 )j7+ 2 (x 2 )Fe//(Xi - X 2 )fj a2 (x 2 )fj ai (x 1 ) , 

-rr , , ^ G+ Q/3 ^ (Xl -X 2 )G /3l/ 3 2;QQ (Xl -X 2 ) 

X4//(xi - x 2 ) = ^ . 

a a z,c a c «i c a 2 

P1P2 

We can see that the second order in V of perturbation theory for S-matrix is equivalent to the first order of perturbation 
theory over effective interaction V e ff- This effective interaction is determined by van der Waals forces (see formula 

<S2»- 

It follows from (|39|l , that the Hamiltonian of interaction of atoms at low-energies is defined by 



V = J dxidx 2 v a i3 ;l s{xi - x 2 )?7+(xi)7)^(x 2 )?7 7 (x2)??5(xi), 



where 

Vafrysix) = i (x 2 (d aS d f } 1 ) - 3(xd Q<5 )(xd,5 7 )) 
This Hamiltonian corresponds to the dipole-dipole interaction of atoms. 



IX. CONCLUSION 



The goal of this paper was to develop a microscopic approach for describing the physical processes in many-particle 
systems in the presence of bound states of particles. To achieve this goal we developed an approximate second 
quantization method for systems with bound states of particles. 

The basic results obtained in this paper are the following: 

1. The Fock space was introduced in the second quantization formalism. In this space the creation and annihilation 
operators of elementary particles x + , \ and their bound states f] + , f] were introduced on an equal status. 

2. The operators of basic physical quantities acting in this space were constructed. These operators include the 
Hamiltonians of interaction between elementary particles and their bound states. 

3. It was shown that in the approximation when the radius of interaction is small the above mentioned Hamiltonians 
transform into the well-known Hamiltonians for Coulomb's and dipole interactions between the particles of various 
kinds. 

4. The non-relativistic quantum electrodynamics of charged and neutral (the bound states) particles was constructed. 

5. The various physical effects including the theory of the van der Waals forces acting between atoms were considered 
as the approbation of the developed formalism. The description of such effects within the usual formalism requires 
more considerable efforts associated with introduction of the interaction between neutral currents of bound states and 
electromagnetic field. 

6. More detailed calculations, which fall out the limits of our approximation should result in appearance of three- and 
many-body interactions between unbound particles and two-body bound states. These many-body interactions should 
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be consistent with pair interaction of original particles. The appearance of these many-body terms is analogous to the 
origin of the terms with photon-photon scattering in transition from the usual standard Lagrangian to the effective 
low-energy Euler-Heisenberg Lagrangian of quantum electrodynamics. 

Especially we would like to emphasize a role of the obtained Hamiltonians, which describe the interaction of quasi- 
neutral particles (the bound states of charged fermions) with electromagnetic field, the elementary particles with 
bound states, and also the bound states with bound states. On the basis of these Hamiltonians one can study 
such phenomena as Bose-Einstein condensation in a gas of excited atoms, the interaction of condensates with an 
electromagnetic field in Bose and Fermi systems. These Hamiltonians can be also the basis of the kinetic theory for 
systems with bound states of particles. 

Finally we would like to stress that the developed method can be easily generalized to the case of bound states 
consisting of more than two particles. The generalization of the offered method for describing the systems with bound 
states of bosons and also of bosons and fermions taking into account the spin of particles can be also performed 
without principal difficulties. 
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For Coulomb's interaction the inequality R » r is equivalent to £ « e because of r = 1 - l —tt, Sq = „'"' : , , n ~ 1 



